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ABSTRACT 

We constrain the linear and quadratic bias parameters from the configuration dependence of the 
three-point correlation function (3PCF) in both redshift and projected space, utilizing measurements 
of spectroscopic galaxies in the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample. We show 
that bright galaxies (M,. < —21.5) are biased tracers of mass, measured at a significance of 4.5cr in 
redshift space and 2.5(t in projected space by using a thorough error analysis in the quasi-linear regime 
(9 — 27 /i~^Mpc). Measurements on a fainter galaxy sample are consistent with an unbiased model. 
We demonstrate that a linear bias model appears sufficient to explain the galaxy-mass bias of our 
samples, although a model using both linear and quadratic terms results in a better fit. In contrast, 
the bias values obtained from the linear model appear in better agreement with the data by inspection 
of the relative bias, and yield implied values of erg that are more consistent with current constraints. 
We investigate the covariance of the 3PCF, which itself is a measurement of galaxy clustering. We 
assess the accuracy of our error estimates by comparing results from mock galaxy catalogs to jackknife 
re-sampling methods. We identify significant differences in the structure of the covariance. However, 
the impact of these discrepancies appears to be mitigated by an eigenmode analysis that can account 
for the noisy, unresolved modes. Our results demonstrate that using this technique is sufficient to 
remove potential systematics even when using less-than-ideal methods to estimate errors. 
Subject headings: large-scale structure of universe - galaxies: statistics - cosmology: observations 



L INTRODUCTION 

Studying the statistical properties of the galaxy dis- 
tribution allows one to probe the structure of overdense 
regions today, learning about galaxy formation and cos- 
mology. We observe significant clumping in this large- 
scale structure (LSS), which is commonly characterized 
by a series of n-point correlation functions (reviewed in 
iPeeblesI I198C0 . Observational evidence is in line with 
predictions of a dar k-energy dominated cold dark mat- 
ter (ACDMl model (iKomatsu et al.ll2009l: Islnchez et al.l 
l2009l : IR,eid eT^MM ^ However, there is a large concep- 
tual hurdle between following the evol ution of mass den- 
sities in gravitational collapse (e.g. iBernardeau et ahl 
l2002f ) and that realized by galaxy positions. A priori, 
there is little reason to believe a one-to-one correspon- 
dence exists between mass overdensities and galaxy posi- 
tions; complex galaxy formation processes such as merg- 
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ing and feedback should have significant contributions. 
For example re cent results froni th e Sloan Digit a ,l Sky 
burvey (SDSS; I York et al.l[200l in lZehavi et all (I200I 
l201Clf ) show clustering varies with galaxy luminosity and 
color. This discrepancy between the observed "light" in 
galaxies relative to the predicted "mass" clustering is of- 
ten described as galaxy-mass bias. 

The parameterization of galaxy-mass bias enables a 
two-pronged approach to probe both cosmology and 
galaxy formation. On one side, we map the cluster- 
ing of galaxies to that of the underlying mass distri- 
bution allowing us to understand and constrain cosmol- 
ogy. Alternatively, the parameterization of the bias itself 
encodes useful information concerning galaxy formation 
processes. This approach distills observational data from 
hundreds of thousands of galaxies available in modern 
surveys, such as the the two-degree field galaxy redshift 
survey (2dFGRS; CoUess et al. 2001) and the SDSS into 
a significantly smaller and more manageable form. 

Most observational evidence exploits the two-point cor- 
relation function (2PCF), the first in the series of n- 
point functions (or equivalently, the power spectrum in 
Fourier space). However, the 2PCF represents only a 
portion of the available information. Measurements of 
higher order moments, such as the three-point correla- 
tion function (3PCF), allow a more complete picture of 
the galaxy distribution. The statistical strength of higher 
order information might rival that of two-point statis- 
tics fSefu satti &: Scoccimarrol 120051 ) . as well as break 
model degenerac ies describing c osmology and ga laxy bias 
(Zhe ng fc Wein berg 2007; Kulka rni et alll2007f ). 

Previous analyses have estimated the 3PCF from 
modern galaxy redshift surveys, inc l uding work on the 
the 2dFGRS fJing & BJmeJ[200l IWang et all [200l 
iGaztahaga et al.. .2005.) and results from SDSS data 
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(iKavo et all [20p: iNichol et all [2001 iKulkarni et alj 
I2007t iGaztafiaga et all 120091: iMarinl I2010D . Related 
higher order statistics have als o been measured for 
these dataset s (Verde ct al. 200"1 iPan fc Szapudil [20051 : 
iHikage et al1[2005t iNishimichi et al.1 12007'^. 

This work is the second of two papers analyzing the 
redu ced 3PCF on SPSS galaxy samples. The first pa- 
per (jMcBride et al.|[2"010[) focused on the details of the 
measurements we analyze here, as well as clustering dif- 
ferences due to galaxy luminosity and color. This pa- 
per utilizes the configuration dependence to constrain 
non-lin ear galaxy-mass bias pa rameters in the local bias 
model ()Frv fc Gaztanagall993D , and the properties of the 
errors necessary for quantitative analyses. 

The local bias model is a simple approach to charac- 
terize galaxy-mass bias. Alternative d escriptions exist 
based on the halo model (reviewed in Coo rav fc ShethI 
l2002f l. which form phenomenological models with a 
wider range of parameters. Two well used formula- 
tions include the halo oc cupation distribution (HOD; 
iBerlind fc Weinberg I2002D and the conditional luminos- 
ity fu nction fCLF: lYang van den Bosch et all 
[200^ . There are formulations for the 3PCF; however, 
the accuracy of the model predictions is not as well de- 
termined as the 2PCF when compared with data (see 
Takada fc Ja"in[ [2001 iWang et all [200l iFosalba et all 
20051) . A significant advantage of a HOD modeling is the 
ability to use well determined measurements of the small 
scales for constraints (the non-linear regime in gravi- 
tational perturbation theory). Understanding the pro- 
jected 3PCF, Qproj, a major component of this work, 
provides a critical link to obtain reliable measurements 
at these smaller scales from observational galaxy sam- 
ples. 

However, by using this simple prescription for galaxy- 
mass bias, we investigate the effects of binning and co- 
variance resolution in a quantitative analysis with a clear 
and simple model where the implications for bias and 
cosmology are better studied for higher order moments. 
An important part of our analysis is comparing results 
from the projected 3PCF with the more commonly used 
redshift space measurements. 

This paper is organized as follows. We discuss the 
SDSS data, simulations, and mock galaxy catalogs in 
|J2| We review the theory and methods of our analy- 
sis in f}3l We constrain the non-linear galaxy mass bias 
parameters in 21 In fJSJ we investigate clustering proper- 
ties contained in the eigenvectors of the 3PCF covariance 
matrix. We perform a detailed examination of the qual- 
ity of error estimation in fJSj We discuss our results and 
compare to related analyses in ijTl Finally, we review our 
main conclusions in fJSj Unless otherwise specified, we as- 
sume aflat ACDM cosmology where = 0.3, VIa = 0.7, 
and Ho — /i 100 kms~^Mpc~^, used to convert redshift 
to physical distances. 

2. DATA 

2.1. SDSS Galaxy Samples 

The SDSS has revolutionized many fields in astronomy, 
obtaining images and spectra covering nearly a quarter 
of the sky by utflizing a dedicated 2.5 meter telescope at 
Apache Point Observatory in New Mexi co (iGunn et al.l 
[T99&. ,2006; .York et al...2000; .Stoughton et al.ll200ir 



Our galaxy samples and details of the measurements 
are fu lly described in a companion paper (jMcBride et al.l 
[2OT0h . Briefly, we use galaxy data with spectroscopi- 
cally det ermined redshifts, d efined as the Main galaxy 
sample (jStrauss et al.l l2002| ). We conduct our anal- 
ysis of clustering measurements usi ng galaxies from 
DR6 (jAdelman-McCarthv et al.ll2008[) . and define sam- 
ples from a refined parent catalog: the New York 
University Va. l ue-Ad dcd Galaxy Catalog (NYU-VAGC; 
iBlanton et al.l l2005[ ). We analyze two samples: a 
BRIGHT sample where Mr < -21.5 and LSTAR 
with —21.5 < Mr < —20.5. We do not analyze 
the FAINT sample p resented in the companion paper 
(jMcBride et al.|[2010D , as the errors suffer from small vol- 
ume effects. We tabulate properties, such as the redshift 
range, number of objects, volume and completeness cor- 
rected number density in Table [T] 

Our absolute r-band magnitudes use the NYU-VAGC 
conventio n defined to represen t values at z — 0.1 (see 
details in IBlanton et al.l 120051) . We note these as Mr 
for simplicity, which refer to Mo.v — 5\ogh. Radial dis- 
tances and absolute magnitudes are calculated using a 
flat ACDM cosmology with = 0.3 and J^a = 1 — ^^m- 

2.2. Hubble Volume Simulation 

To estimate the clustering of mass in the late time 
ACDM cosmology, we analyze cosmological A^-body sim- 
ulations. We use the Hubble Volume ( HV) simulations 
(jColberg et al.ll2Q00l: lEvrard et aL[|200a) that were com- 
pleted by the the Virgo Consortium. We utilize the light- 
cone output with ACDM cosmology: (fim = 0.3, ^Ia — 
0.7, h = 0.7, 0-8 = 0.9), where Ho = hlQQ kms'^Mpc'^ 
The HV simulation consists of 1000'^ particles in a box 
of (3000 h~^Mpc)^ volume, resulting in a particle mass 
of rripart — 2.2 x 10^^ h~^MQ. The particles start from 
an initial redshift of Zinu — 35, and are evolved to the 
current time using a Plummer softened gravitational po- 
tential with a softening length of 0.1 /i~^Mpc. 

We use the sam e simulation output a s presented in our 
companion paper (jMcBride et al.ll2010D . Here, we briefly 
review our postprocessing of the simulation data for com- 
pleteness. We include redshift distortions in the mass 
field by distorting the position according to the peculiar 
velocity of the dark matter particle. We trim particles to 
match the identical volume of the corresponding SDSS 
samples, including the non-trivial angular geometry of 
SDSS data. Finally, we randomly downsample the num- 
ber of dark matter particles to make the computational 
tin ie of the analysis mo re manageable (discussed further 
in IMcBride etaHIMol) . 

2.3. Mock Galaxy Catalogs 

We analyze mock galaxy catalogs created to match 
some of the SDSS galaxy data. These mock catalogs 
were constructed from 49 independent A'^-body simula- 
tions, initiated with different random phases and evolved 
from a single cosmology: (ilm — 0.27, n\ = 0.73, h — 
0.72, CTs — 0.9). While these differ slightly from our as- 
sumed cosmology for the data, we expect the differences 
to be minor, with no significant implications on our anal- 
ysis. Each of the 49 realizations had randomized phases 
where the initial conditions were generated from 2nd or- 
der Lagrangian perturbation theory ((Scoccimarro. .19981 : 
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Sample 


Specifics of SDSS galaxy samples 

Absolute J , .r. Volume Number of 
Magnitude h-^Gpc^ Galaxies 


Density 
10-3 h^Mpc-^ 


BRIGHT 
LSTAR 


Mr < -21.5 0.010 - 0.210 
-21.5 < Mr < -20.5 0.053 - 0.138 


0.1390 37,875 
0.0391 106,823 


0.272 
2.732 



Table 1 

The magnitude range, redshift limits, volume, total number of galaxies, and completeness corrected number density are shown for the 
galaxy samples constructed from the SDSS DR6 spectroscopic catalog. We selected these samples by cuts in re dshift and corrected 
(K-correction and passive evolution) absolute r-band magnitude to create volume-limited selections. See details in lMcBride et aLl 1120101 ). 



iCrocce et al.ll200^ . These simulations each c onsist of 
640^ particles that we evolved using Gadget2 (|Springell 
l2005f ) from an initial redshift of Zi = 49 to the present 
epoch. The box side- length of 1280 h~^Mpc contained 
enough volume to exactly match the brightest galaxy 
sample after applying the SDSS geometry. These sim- 
ulations have been used in various o ther studies (e.g. 
[Tinker et al1[200l iMamsra et al.|[2QTol ). 

The galaxy mocks were created by populating dark 
m atter halos wi t h gala xies by applying the HOD model 
in [Tinker et alj ()2005[ ) with parameter values defined 
to represent Mr < —21.5 and as — 0.9. The ha- 
los were id entified using a friends-of-friends algorithm 
(j Davis et a l. 1985) applying a linking length of 6 = 0.2 
in units of the mean interparticle separation. The least 
massive halos contained 33 particles, a minimum mass 
capable of hosting the faintest galaxies in the BRIGHT 
galaxy sample. Given the mass resolution of these sim- 
ulations, less massive halos necessary to host galaxies in 
the LSTAR galaxy sample could not be identified. There- 
fore, we could only obtain reliable mock galaxy catalogs 
corresponding to the BRIGHT sample. 

3. THEORY & METHODS 

The n-point correlation functions remain a standard 
descri ption of the co mplexity seen in large-scale structure 
(LSS: IPeebleslll980l ). In terms of the fractional overden- 
sity (S) about the mean density (p), 



p{x) 

p 



1 



(1) 



we characterize the two-point correlation function 
(2PCF) and three-point correlation function (3PCF) as: 

ari2) = (<5(fi)<5(x2)) . (2) 

C(''l2,''23,''3l) = {S{X1)6{X2)S{X3)) . (3) 

We make the standard assumption of a homogeneous and 
isotropic distribution, and report clustering amplitudes 
dependent on the magnitude of the separation vector, 

e.g. ri2 = |xi -fzl- 

Motivated by the hierarchical ansa tz (iPeeblesI 119801) 
and gravitational perturbation theory ([Bernardeau et aLl 
[200 a we use the reduced 3PCF: 



Q{ri2,r23,r3i) 



C{ri2,r23,r3i) 



'^12f23 + $23^31 + ^31^12 



(4) 



This "ratio statistic" remains close to unity at all scales, 
and to leading order is insen sitive to both time evolu tion 
and cosmology (reviewed in [Bernardeau et aL|[2002D . 

Redshift distortions impact measurements of cluster- 
ing by altering the line-of-sight radial distant estimate, 



as we are unable to distinguish the gal axy's peculiar ve - 
locity from the Hubble flow (reviewed in [Hamilton[[1998l ). 
We refer to the theoretical non-distorted distances as real 
space, commonly denoted with r. Distances that include 
the redshift distortion (e.g. observational distances) are 
in redshift space, denoted with s. We decompose the red- 
shift space distance into line-of-sight (tt) and projected 
separation (rp) such that s = (tt^ -t-r^)^/^. With this sep- 
aration, the anisotropic distortion is primarily contained 
in the tt coordinate. 

We minimize the impact of redshift distortions by es- 
timating the correlation function binned in both rp and 
TT and integrate along the line-of-sight resulting in the 
projected correlation function (Davis & Peebles., 19831) : 



Wp{rp) = 2 ^(rj,,7r)d7r 



(5) 

The projected 3PCF and its reduced form have analogous 
definitions: 

Cproj {rpl2 , f'p23 ,rp3l) = 

C(^pl2,7'p23,?'p31,7I"l2,7r23)d7ri2d7r23 , (6) 



]proj {rpl2, rp23 , Tp31 ) 



Cproj (?'pl2 , rp23 , r-p3i ) 



Wpl2Wp23 + Wp23Wp3i + Wp3iWpi2 



(7) 



The measurements we analyze set TTmax = 20 /i^^Mpc. 
We find this sufficiently deep to recover correlated struc- 
ture to minimize redshift distortions, but not overly ex- 
pe nsive to calculate (see detailed discussion in appendix 
of IMcBri de et aLl[2010l) . 

The full 3PCF is a function of three variables that 
characterize both the size and shape of triplets. We pa- 
rameterize the 3PCF by (ri, r2, 9), where ri and r2 rep- 
resent two sides of a triangle (simplified notation from 
ri2 and r23), and 9 defines the opening angle between 
these sides. However, our measurements are estimated 
in bins defined by (^12, r-23, rsi). We convert r.^i to 9 us - 
ing the cosine rule (as detailed in [McBride et aLl[2010D . 
The 3PCF remains sensitive to the exact choice of bin- 
ning sch£me;_wliichrajijM^ 

nal ([Gaztafiaga &: Scoccimarrol [2005[ : [Marin et all [2008[ : 
[McBride et al.[[2010D . We choose a bin-width as a frac- 
tion, /, of the measured scale, r, such that = / x r 
and a bin at r represents (r — r + ■^). 

We always use the reduced 3PCF as a function of three 
variables, Q{ri,r2,9), but we simply our notation by 
sometimes referring to it as Q{9) or even Q. If the ampli- 
tude of Q{9) varies significantly with 9, we refer to this 
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as strong configuration dependence, in contrast to little 
or no variation for a weak configuration dependence. We 
define the scale of triangles by ri , and choose configura- 
tions such that r2 = 2ri. This results in rg, varying in 
size from ~ r2 ~ vi when 6 ~ Q to = r2 + ri when 
61 = TT. 

3.1. Galaxy-Mass Bias 

We can consider galaxies to be a biased realization 
of the ACDM mass fi eld. In the local bias model 
(|Frv fc Gaztanagall993D . the galaxy over-density, 5g, can 
be connected to the mass over-density, dm, by a non- 
linear Taylor series expansion: 



E 



k\ " 



h ,2 



(8) 



This relation describes the mapping between galaxy and 
mass by simple scalar values, to second order: the linear 
(hi) and quadratic (62) bias. 

With measurements on galaxy n-point correlation 
functions, the clustering of galaxies is linked to mass clus- 
tering via the bias parameters. The 2PCF can be used 
to constrain the linear bias by equating the correlation 
function between galaxies, ^g, to that of dark matter, 
such that 

Ur)^blUr). (9) 

The 3PCF is the lowest order correlation function that 
shows leading order sensitivity to the quadratic bias 
term. The analog to © for the connected 3PCF is writ- 
ten 



Cg{ri2,r23,r3i) = blCm{ri2,r23,r3i) + 

b\h2 [C12C23 + Cl2^31 + ^31?2 



(10) 



where ^12 = 6n(''i2), etc. This simplifies for the reduced 
3PCF where we denote the bias parameters ds B = hi 
and C &2/&1: 



1 



Qg{ri2,r23,0) = -[Q™(ri2,r23,^)+C] 



(11) 



We have changed notation slightly in (ITT]) , replacing r^i 
with 0, the opening angle between the two sides ri2 and 
r23, as we discussed above. 

A multiplicative factor such as B can dampen {B > 1) 
or enhance {B < 1) the configuration dependence of Qm 
as seen from the galaxy distribution, whereas the value 
of C will produce an offset. We see that B and C are par- 
tially degenerate in this model. If Q{9) shows no config- 
uration dependence, two parameters are used to describe 
a shift in amplitude. However, this degeneracy can b e re- 
move d when the 3PCF exhibits a shape dependence (jFrvl 
119941 ). Even with the degeneracy broken, the values of 
B and C could show a strong correlation. 



for each realization (k) and bin (i) given a mean value 
(Qi) and variance (cr|) for each bin over all realizations. 
We use Q as a general placeholder for any measured 
statistic (2PCF, 3PCF, etc). 

We construct the normalized covariance matrix using 
the standard unbiased estimator: 



C 



N 



1 ^ 



(13) 



k=l 



Equation [13] assumes that each realization is indepen- 
dent. In practice, a number of mock galaxy catalogs can 
be used to make this a tractable approach. If mock cata- 
logs appropriate to the galaxy sample are not available, a 
covariance matrix can be estimated from the data itself, 
such as the commo nly employed jackknife re-sampling 
(|Lupton et al.ll2001[ ). Since jackknife samples are not in- 
dependent realizations, we compute the covariance by: 



C. 



Uack) _iN~ If 



N 



N 1 ^ 
C - - - i Va'^A'^ 



(14) 



k=l 



where Cy denotes the typical unbiased estimator of the 
covariance when computed on N jackknife samples. 

3.3. Eigenmode Analysis 

We constrain galaxy-mass bias parameters using the 
information in the full covariance matrix. We utilize 
an eigenmode analysis (Scoccimarro 2000), an equivalent 
method to a principal component analysis (PGA) on the 
measurement covariance matrix. This method was tested 
in detail for th e galaxy-mass bias of the 3PCF us ing sim- 
ulated data in iGaztanaga fc Scoccimarrol (j2005| ) . 

The basic idea is to isolate the primary contributing 
cigcnmodcs of the reduced 3PCF based on the structure 
of the normalized covariance matrix. This allows one 
to trim unresolved modes and perform a fit in a basis 
which minimizes the non-Gaussianity of the residuals. To 
summarize, the covariance matrix can be cast in terms 
of a singular value decomposition (SVD), 

- C/ E ; = X^6,, . (15) 

where Sij is the Kronecker delta function making S a 
diagonal matrix containing the singular values, Xf. The 
matrices U and V are orthogonal rotations to diagonalize 
the covariance into E where V'^ denotes the transpose of 
V. 

Applying the SVD to the covariance matrix yields a 
rotation into a basis where the eigenmodes are uncorre- 
lated (i.e. the covariance matrix becomes diagonal) . The 
resulting rotation matrix can be directly applied to our 
signal forming the Q -eigenmodes, 



(16) 



3.2. Estimating the Covariance Matrix 

We measure the correlation between measurements by 
empirically calculating the covariance matrix. Given a 
number of realizations, N, a fractional error on Q can be 
written as 



The singular values provide a weight on the importance 
of each eigenvector. Specifically, a multiplicative factor 
of 1/Xf is applied when C is inverted. With this feature 
in mind, we define the signal-to-noise ratio as 



(17) 
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We note that this S/N estimate is a lower bound on the 
true S/N due to the SVD. To remove noise and avoid 
numerical instabihties, we tri m eigenmodes correspond- 
ing to low singular values. iGaztafiaga fc Scoccimarrol 
([2005') suggest keeping eigenmodes resolved better than 
the sampling error in the covariance matrix. Since our 
covariance matrices are normalized (i.e. the diagonal el- 
ements are set to one), the singular values are directly 
related to sampling error, and we require the so-ca lled 
"dominant modes" (jGaztanaga fc Scoccimarrol 120051 ) to 
satisfy: 

\^>V2/N, (18) 

where N refers to the number of samples used to estimate 
the covariance matrix. 

The advantage to using this eigenmode analysis for 
fitting is threefold. First, it correctly incorporates the 
correlation between measurement bins. Second, by per- 
forming the fit in the rotated basis of the eigenmodes, the 
residuals of the fit are more Gaussian and the degrees of 
freedom are properly addressed (e.g. 3 eigenmodes really 
only fits over 3 numbers). Finally, using only dominant 
modes removes artifacts due to noise in the estimated 
covariance matrices. For example, when using the full 
covariance but not trimming any modes, noise can cause 
a fit to converge on incorrect values with artificially small 
errors (and falsely high S/N) . This effect becomes worse 
as the covariance becomes less resolved. Conversely, fit- 
ting over dominant eigenmodes helps to eliminate any 
roblems from unresolved parts of the error estimation 
Gaztanaga &: Scoccimarrol l2005. Figure 13), and has the 
benefit of dealing with singular covariance matrices. 

4. GALAXY-MASS BIAS 

We want to constrain the galaxy-mass bias described 
by ([8]) using the full configuration dependence of the 
reduced 3PCF in the quasi-linear regime. For the 
galaxy data, we use measurements in both redshift space, 
Qz(d)i and projecte d space, Qproj(0) as presented in 
iMcBride et al.l (|2010D . We estimate the bias parameters 
by comparing to mass estimates obtained from dark mat- 
ter particles in the HV simulation (see §ss:hv). We expect 
redshift distortio ns to affect the bias relat ion, which we 
partially neglect (jScoccimarro et afl 11999( 1 . In particu- 
lar, we account for the effects of redshift distortions by 
applying a distortion distance to the dark matter parti- 
cles based on their velocities for our mass measurement. 
However, this is not completely sufficient as redshift dis- 
tortions alter the bias re l ation in (fTT|) . especially for Qz 
(jScoccimarro et aLlll999L l2001h . We expect that Qproj 
will be predominantly unaffected and roughly equivalent 
to r eal space mea surements for this parameterization (see 
e.g. lZheng|[20Ql . 

We restrict our analysis to scales above 6 /i~^Mpc, 
corresponding to Q{ri,r2,0) with ri = 6 and 9 h~^Mpc 
with configurations having r2 = 2ri. W e let 9 vary be- 
tween and TT using 15 bins, as detailed in lMcBride et al.l 
(|2010( ). We investigate galaxy-mass bias in two sam- 
ples where the covariance is well determined: BRIGHT 
(Mr < -21.5) and LSTAR (-21.5 < Mr < -20.5) as 
listed in Table [TJ We remove the least significant eigen- 
modes during the fit by applying the criteria in ([T5)) . 
We discuss possible effects of using a different number 
of modes in ^ For each galaxy sample, we perform 
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Figure 1. Constraints on the galaxy-mass bias parameters using 
the Mr < —21.5 galaxy sample and the HV simulation for mass 
estimates. The left column corresponds to fits using Qz{9) (redshift 
space) with the right column fit using Qproj {S) (projected space). 
The top and bottom panels represent individual fits with triangles 
of ri = 6 and 9 ?t"^Mpc as indicated. The middle panels are a 
joint fit using both triangles. There are two points of comparison 
marked: an unbiased result with {B = 1.0; C = 0.0) and only non- 
linear bias (B = 1.0; C = —0.3). The contours denote the 1,2 and 
3(T levels from the Ax^ distribution of two parameters. 

six independent fits: a series of three different scales 
for measurements in both redshift and projected space. 
Wc use the full configuration dependence for triangles 
with ri = 6 and 9 /i^^Mpc, as well as a joint fit using 
both scales. For the joint fit, we estimate the full com- 
bined covariance matrix to correctly account for overlap 
and correlation and use the same eigenmode analysis. 
This changes the number of available modes from 15 in 
the individual fits to 30 modes for the combined joint 
fit. We estimate the covariance for these samples us- 
ing 30 jackknife samples, where our our jackknife regions 
have equal unmasked area on the sky and use the full 
redshift distribution of the observational galaxy sample 
(|McBride et al.llMol) . 

4.1. Constraining Non-Linear Local Bias 

We constrain the galaxy-mass bias using a maximum 
likelihood approach by calculating a simple statistic 
where the likelihood C cx exp(— x'^/2) and 

A,^ ^' ^' . (19) 
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Figure 3. The reduced 3PCF for the Mr < —21.5 sample show- 
ing the mass scaled to the "best fit" galaxy-mass bias parameters. 
The top two panels correspond to redshift space, and the bottom 
two to projected space. From left to right, the scale of the triangle 
increases as noted. The red (dashed) line represents an individual 
fit only to that triangle scale, and the blue (dotted) line shows a 
joint fit between both scales. 



Figure 2. Analogous to Figure [T] but for -21.5 < Mr < -20.5 
galaxies. 

We determine the theoretical model, Q^*\ by scaling 
the mass measurement from the HV simulation, Qm, 
with bias parameters B and C as per (ITT]) . We evalu- 
ate £ on a grid using the ranges: B = 0.1 . . .3.0 and 
C = —1.5 ... 1.5 with a step-size of 0.01. We tested for 
discrepancies using a factor of 10 finer spacing between 
grid elements with no significant differences to the fitted 
results. 

We first examine the BRIGHT sample {Mr < -21.5), 
with the likelihood space of the six 2-parameter fits dis- 
played in Figure [TJ We include contours for Gaussian 

I, 2 and 3cr levels which identify regions of probabilities 
for 68.3,95.5 and 99.7%. We calculate these from the 
Ax^ distribution for a 2-parameter fit (i.e. two degrees- 
of-freedom), with corresponding values of 2.3, 6.2, and 

II. 8 from the best fit value. We include two refer- 
ence points for comparison, the unbiased result where 
{B = 1.0; C = 0.0) along with a potential negative 
quadratic bias term accounting for the entire galaxy bias 
(B = 1.0; C = ~0-3) ( similar comparison to Figure 5 in 
iGaztafiaga et al.ll2005l) . 

We can clearly see the degeneracy between B and C 
in Figure [TJ visible as the elongated diagonal contour. 
Larger values of B remain likely with larger values of 
C, consistent with our expectation of degeneracy by in- 
specting the bias relation in (ITT]) . The size of the errors 
are notably larger for projected space measurements, as 
well as lower values for the overall S/N. This results 
from the larger uncertainties in the projected measure- 
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Figure 4. Like Figure|3]but for the -21.5 < Mr < -20.5 sam- 
ple. The reduced 3PCF showing the mass scaled to the "best fit" 
galaxy-mass bias parameters. The top two panels correspond to 
redshift space, and the bottom two to projected space. From left to 
right, the scale of the triangle increases as noted. The red (dashed) 
line represents an individual fit only to that triangle scale, and the 
blue (dotted) line shows a joint fit between both scales. 
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ments (IMcBride et aUlMoh . Since the scale Vp repre- 
sents a projection that incorporates larger scales (deter- 
mined by the line-of-sight integration iTmax), projected 
measurements are more sensitive to the dominant uncer- 
tainty from cosmic variance that increases with scale. In 
all cases, the unbiased {B — l;C — 0) model is excluded 
at greater than a 2a level. To see the success of the fit 
"by eye" , we plot the 3PCF for dark matter, galaxies and 
best fit scaled model for this sample in Figure [3] Both 
the "individual" fits and "combined" joint fit produce 
models that well match the data. 

Next, we fit the galaxy- mass bias parameters using the 
LSTAR sample (-21.5 < Mr < -20.5). This sample 
spans a unit bin in magnitude, and consists only of galax- 
ies fainter than the previous bright sample. The results 
of the fit with likelihood contours arc shown in Figure O 
The uncertainties appear reduced in size - a striking dif- 
ference with respect to the BRIGHT sample in Figure [TJ 
In addition, the slope of the "line of degeneracy" between 
B and C has shifted. We reason that this is in part due 
to the increased statistical significance of the larger sam- 
ple, as both the measurements and covariance are better 
resolved. Due to the higher number density of galaxies, 
we re-measured the 3PCF using a finer binning scheme 
(fractional bin-width oi f — O.I as opposed to / = 0.25, 
see comparison in Appendix With the finer binning, 
we see a stronger configuration dependence, which will 
alter the degeneracy between B and C. We note that 
many of the best fit B values appear smalle r, which w e 
expect for a fainter sample (Zehavi et al.. ,20051 I2Q10D . 
The same line of reasoning suggests that the "unbiased" 
model {B = 1]C = 0) should be more likely to fit. 

As before, we plot the respective best fit model in com- 
parison with the dark matter and galaxy 3PCF in Fig- 
ure |4l There is a smaller difference between HV (mass) 
and galaxy measurements, as this sample is fainter. We 
notice some noise of the HV measurement for Qproj , mak- 
ing the model not quite as smooth. We note that by eye, 
Qz on larger scales indicates a slight bias for the com- 
bined fit, with the model undershooting the data and la 
uncertainties. Significant off-diagonal structure in the 
covariance matrix can produce a fit where "chi-by-eye" 
suggests a poor fit. Since the ri = 6 h^^Mpc mea- 
surements in Qz have much smaller errors, these scales 
drive the fit making measurements with ri = 9 /i~^Mpc 
appear a poor match to the "best fit" model. 

We summarize the results of our two parameter con- 
straints for the BRIGHT and LSTAR sample in Tabled 
The BRIGHT sample (Mr < -21.5) represents galax- 
ies with r-band magnit udes significantly bri ghter than 
where Mr 20.4 (jBlanton et al.ll2003[ ). We typ- 
ically consider L, galaxies to have a linear bias, i.e. 
where B 1, and we might expect this brighter sam- 
ple to have a larger B value. The constraints from 
projected measurements appear to follow this logic; the 
best fit values on Qproj in the fainter LSTAR sample 
(-21.5 < Mr < -20.5) are lower with B ~ 1. Redshift 
space measurements, Qz, appear consistent with B ^ 1 
for all fits, but at the same time values of C are lower, 
reflecting the degeneracy of the B and C parameters. 
The reduced xt values show an acceptable fit in almost 
all cases; the exceptions are the two Qz fits using the 
ri = 6 h~^Mpc triangles for the LSTAR sample. Con- 
sequently, the joint fit appears to be the poorest match in 



Figure m The Ax^ in Table [5] displays the likelihood an 
unbiased model is from the best fit parameters. We find 
an unbiased model is ruled out for the BRIGHT galaxy 
sample at greater than A.Sa in redshift space and 2.6a in 
projected space. We cannot conclude the same for the 
LSTAR sample, which is largely consistent with an un- 
biased model. We generally consider bri ght galaxies to 
be more biased (jZehavi et all 120051 . [2010f ). The LSTAR 
sample is a magnitude bin around and fainter than 
the BRIGHT sample, and we expect a better consistency 
with the unbiased model. 

4.2. Non-zero Quadratic Bias? 

With our two parameter likelihood space, we can inves- 
tigate the statistical significance of a non-zero quadratic 
bias term (62 which is encapsulated in C = 62/^1)- We 
use the same configuration dependence of the 3PCF, and 
the measured covariance, but restrict the two parameter 
fit such that C = 0. We evaluate the best fit B, the qual- 
ity of the fit (via the reduced xt) as well as the Ax^ for 
the best two-parameter fit, which we present in Table [3l 
For the BRIGHT sample, we notice the B values are 
equivalent across both Qproj and Qz for the same scales 
on the same sample. Since we removed the degeneracy 
(as C is zero), this behavior makes sense and in agree- 
ment with the measurements. We note that the typical 
B values are larger for the BRIGHT sample, and lower 
for the fainter LSTAR sample. For Qproj our constraints 
find little statistical significance for a non-zero quadratic 
bias term; the likelihood difference is small and a lin- 
ear bias term is sufficient to quantify the bias for both 
the BRIGHT and LSTAR samples. Overall, measure- 
ments in redshift space (Qz) more strongly suggest that 
C ^ 0, especially when using the smaller scale triangles 
(n = 6 h-^Mpc). 

4.3. Relative Bias 

The relative bias characterizes the relative clustering 
strength between different galaxy samples - an alterna- 
tive to the "absolute" galaxy-mass bias constrained pre- 
viously. Relative bias is insensitive to cosmology and 
does not require assumptions to determine mass clus- 
tering. We can use the relative bias to check consis- 
tency with linear and quadratic bias parameters obtained 
above in ^ For the 2PCF, the relative bias is simply: 

€1 = , (20) 

where ^ can refer to redshift or projected space measure- 
ments. 

We show 6|.7; from the 2PCF in Fi gure m using 
the linear bias parameters obtained from the best two- 
parameter fit (i.e. Table Both redshift space and 
projected measurements agree and produce a flat rela- 
tive bias, even at non-linear scales below a few ft,~^Mpc. 
Two obvious discrepancies arise when comparing obser- 
vational data to "best flt" values. First, neither redshift 
nor projected space fits appear to match data. Earlier 
we noted a substantial degeneracy between the linear 
and quadratic bias terms. The quadratic bias term is 
accounting for more of the clustering bias when we con- 
strain with Q{9), which isn't noticeable in the 2PCF. 
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Measurement Scales ( h -"^Mpc) 



Galaxy-Mass Bias Parameters from SDSS 



B 



C 



D.o.F. unbiased Ax 



BRIGHT-z 
BRIGHT-proj 
BRIGHT-z 
BRIGHT-proj 
BRIGHT-z 
BRIGHT-proj 



LSTAR-z 
LSTAR-proj 
LSTAR-z 
LSTAR-proj 
LSTAR-z 
LSTAR-proj 



6-18 
6-18 
6-27 
6-27 
9-27 
9-27 



+U.ii 



0. 
+0. 

0. 
+0. 

+0. 

0. 
+0. 

0. 
+0. 
0. 

101°; 

+0. 

0. 
-1-0. 

0. 
-1-0. 
0. 

031°; 



-0.25l^-"« 
-0.03t° 
-0.24« 
-0.06t° 

-)+0 
-0 

-0.02t° 



-0.22"* 



1.48 


6-2 


118.86 (10.7(t) 


0.78 


6-2 


16.43 (3.6(t) 


0.83 


9-2 


132.54 (11.3tT) 


0.45 


10-2 


18.70 (3.9<t) 


0.60 


4-2 


26.90 (4.8cr) 


0.34 


5-2 


9.44 (2.60-) 


13.47 


3-2 


28.07 (4.9cr) 


0.85 


4-2 


3.00 (1.2(7) 


3.22 


5-2 


45.85 (6.5cr) 


1.07 


7-2 


5.86 (1.9(7) 


0.07 


3-2 


2.37 (1.0(7) 


1.75 


4-2 


1.93 (0.9(7) 



6-18 
6-18 
6-27 
6-27 
9-27 
9-27 



96: 
03j 
04"* 



-0.22 
-0.01 
-0.30 
-0.14 
-0.07 
-0.09 



Table 2 

The two-parameter best fit galaxy-mass bias parameters, using l lllll with the configuration dependence in the reduced 3PCF from SDSS 
DR6 galaxy samples in comparison with dark matter clustering from the Hubble volume simulation. The fits are performed separately on 
two galaxy samples BRIGHT {Mr < —21.5) and LSTAR, (—21.5 < Mr < —20.5) using measurements in redshift space (denoted with "z") 
as well as projected space ("proj"). The second column lists the range of scales used for the respective fit. The errors are marginalized la 
bounds calculated by the range within Ax^ < 1 from the best fit value. The quality of the best fit value is stated with the reduced 
chi-square = X /D.o.F. The degrees of freedom (D.o.F.) correspond to the number of eigenmodes used minus the number of 
parameters (2forallthesefits). The last column lists the Ax'^ value to quantify the likelihood of an "unbiased" model matching the data, 
i.e. (B = 1;C = 0), with a likelihood expressed in the number of a from by the standard Gaussian assumption for the Ax^ distribution. 



Galaxy-Mass Bias without Quadratic Term 



Measurement 


Scales ( /i^^Mpc) 


B 




D.o.F. 


Ax'^ from best fit 


BRIGHT-z 


6-18 


1.34 


+U.U4 
-0.04 


2.42 


6-1 


6.16 


(2.0(7) 


BRIGHT-proj 


6-18 


1.31 


+0.11 
-0.09 


0.63 


6-1 


0.04 


(0.0(7) 


BRIGHT-z 


6-27 


1.30 


+0.03 
-0.03 


2.23 


9-1 


12.01 


(3.0(7) 


BRIGHT-proj 


6-27 


1.27 


+0.09 
-0.07 


0.42 


10-1 


0.16 


(0.1(7) 


BRIGHT-z 


9-27 


1.24 


+0.06 
-0.06 


1.85 


4-1 


4.35 


(1.6(7) 


BRIGHT-proj 


9-27 


1.25 


+0.11 
-0.09 


0.26 


5-1 


0.01 


(O.Ocr) 


LSTAR-z 


6-18 


1.21 


+U.0!) 
-0.04 


8.68 


3-1 


3.90 


(1.5(7) 


LSTAR-proj 


6-18 


1.11 


+0.07 
-0.06 


0.57 


4-1 


0.01 


(0.0(7) 


LSTAR-z 


6-27 


1.23 


+0.04 
-0.04 


4.59 


5-1 


8.70 


(2.5(7) 


LSTAR-proj 


6-27 


1.15 


+0.07 
-0.07 


1.02 


7-1 


0.76 


(0.4(7) 


LSTAR-z 


9-27 


1.08 


+0.06 
-0.05 


0.14 


3-1 


0.21 


(O.lcr) 


LSTAR-proj 


9-27 


1.11 


+0.09 
-0.08 


1.25 


4-1 


0.25 


(0.1(7) 


Table 3 



The single-parameter best fits for galaxy-mass bias using {TTJ where wo constrain C = 0. We fit the configuration dependence of the 
reduced 3PCF from SDSS DR6 galaxy samples in comparison with dark matter clustering from the Hubble volume simulation. Fits are 
performed separately on two galaxy samples BRIGHT {Mr < —21.5) and LSTAR (—21.5 < Mr < —20.5) using measurements in redshift 
space (denoted with "z") as well as projected space ("proj"). The second column lists the range of scales used for the respective fit. The 

errors are marginalized Icr bounds calculated by the range within Ax^ < 1 from the best fit value. The quality of the best fit value is 
stated with the reduced chi-square xi = X^ /D.o.F. The degrees of freedom (D.o.F.) correspond to the number of eigenmodes used minus 
the number of parameters (in this case, just one). The last column lists the Ax^ value to quantify the difference in likelihood of this 
model with C = compared with the best fit of a two-parameter fit (i.e. Table [2)|- 



This suggests we underpredict values of linear bias, ei- 
ther just for the BRIGHT sample or in unequal portions 
for both. Second, there is a significant difference be- 
tween these two estimates given the same galaxy sam- 
ples, although the projected measurement appears closer 
to agreement. 

Let us consider the relative bias of the reduced 3PCF. 
Since Q is proportional to 1/B, we define the relation 

,(3) ^ QLSTAn_ (21) 
WBRIGHT 



figuration dependence measurements of Q{9) used for 
constraining B and C . Looking at Figure [51 we see an 
obvious difference with respect to the 2PCF: the much 

(3) 

larger uncertainties. However, the predicted 6^^; from 
the galaxy-mass bias constraints appear much more con- 
sistent with the measurements, as opposed to the 2PCF. 
The 3PCF results agree with the observational data, and 
show a much smaller discrepancy between redshift and 
projected space. The quadratic bias term (C) can prop- 
erly account for the clustering difference that was missing 
in the relative bias of the 2PCF. 



Figureinipresents the relative bias of our DR6 galax ies for 
the equilateral 3PCF {Qeq) (jMcBride et aLl^Olol ). We 
note that Q^q is related but not identical to the con- 



4.4. Implications for Cosmology: (Tg 
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Figure 5. The relative bias 6*.^', = ^/(bright /^LSTAR using 
measurements in redshift (red 'x' symbols) and projected space 
(blue diamonds). We calculate the uncertainties by propagating Icr 
values from the 2PCF. The dotted and dashed lines display results 
from the best fit bias terms at the largest scales (9 — 27 h~^Mpc). 
The bold lines indicate values from the two-parameter fit (Table[2|, 
and the faint lines show the best linear fit (Tabled. 



1 








> 




' ' 1 






























"rom 












from 










Redshift Best Fit 




1 


Projected Best Fit 

1 



12 5 10 



r (Mpc / h) 

Figure 6. Analogous to Figure[5]but for the 3PCF. The relative 
bias b^Ji = Qlstar/Q BRIGHT using measurements of equilateral 
3PCF in redshift (red 'x' symbols) and projected space (blue dia- 
monds). We calculate the uncertainties by propagating Icr values 
from the 3PCF. The dotted and dashed lines display results using 
the best fit bias terms at the largest scales (9 — 27 /i~^Mpc). The 
bold lines indicate values from the two-parameter fit (Table[2]l, and 
the faint lines show the best linear fit (quadratic bias constrained 
to be zero; Table |3j. 

Better understanding galaxy-mass bias, or at the least 
accurately parameterizing it, allows one to "calibrate 
out" the effects of galaxies and infer properties of the 
underlying mass distribution to constrain cosmology. We 
can use our estimates of bias to probe the mass variance 
in spheres of 8 /i~^Mpc radius, a common normalization 
of the amplitude of the matter power spectrum, P(k). 



The theoretical erg is linearly extrapolated from a very 
early epoch until today, 

k'^dk 

al=4n W^{k,R = 8 h'^Mpc)PH^{k)—-. , 

(22) 

where W{k,R) is a top-hat window function in Fourier 
space for mode k and smoothing radius R and Pun (k) is 
the linear power spectrum. 

In terms of our fitting formula on the 3PCF in ([TT|) . 
we expand the bias relation for the 2PCF to highlight its 
dependence on erg 

^.M = S'(^)'UW- (23) 

Formally, the mass 2PCF already encodes a value of erg. 
As scales linearly with a change in the square of erg , we 
include an explicit scaling factor to account for a differ- 
ence in CTg between the underlying mass of the observed 
galaxy distribution and that assumed in our estimate of 
mass clustering from A^-body results. In our case, we use 
dark matter from the HV simulation where erg = 0.9, ex- 
plaining the denominator on the right hand side of ()23p . 
We can see that an incorrect assumption of erg in the 
estimate of mass will directly translate into a different 
value of the best B describing galaxies. Even if we use 
the above relation, (|23l) , B and erg are completely degen- 
erate when solely considering the 2PCF. 

By using the additional information available in the 
configuration dependence of the reduced 3PCF, we ob- 
tain a value of B that is independent of erg , and breaks 
the degeneracy between the two parameters. Formally, 
this is only true to leading order, as loop corrections in 
Q{0) will add cosmological dependence which we neglect 
in this analysis. 

With an independent value of B from (fTT|) . we esti- 
mate erg by utilizing the 2PCF in (f23l) . Ideally, we could 
construct a three-parameter fit to jointly constrain B, C 
and erg (e.g. iPan fc Szapudi|[2005l on 2dFGRS data). Or 
as an further extension, we could jointly fit over several 
samples, since they each have the same underlying erg. 
However, this additional complexity is beyond the scope 
of this analysis as our current uncertainties would yield 
poor constraints on erg. We simply estimate the value of 
erg implied by best fit bias parameters. We restrict this 
estimate to the largest scale triangles (ri = 9 h~^Mpc) 
to ensure we approach the linear regime (i.e. the scales 



Implied values of erg 


Measurement 


Scales ( ft "1 Mpc) 


B 


eg 


BRIGHT-z 
BRIGHT-proj 


9-27 
9-27 


1 24+"-"« 
-^■^^-0.06 

1 25+"-" 
-^■^°-0.09 


0.96-1.13 
1.02-1.12 


LSTAR-z 


9-27 




0.88-0.97 


LSTAR-proj 


9-27 


1 -,-,-1-0.09 

-'-■-'--'--0.08 


0.83-0.97 



Table 4 

We use galaxy-mass bias constraints from the configuration 
dependence of the 3PCF, Q{6), with me asur ements of the 2PCF, 
to estimate the implied values of erg via II23I I. We use the largest 
triangle configurations for our two samples, and the 1-parameters 

constraints on B. The range of erg does not represent formal 
uncertainties; we calculate values from the range of uncertainties 
stated in B, neglecting additional errors from the 2PCF. For 
reference, WMAP- 5 (with SN and B AO ) suggest erg = 0.82 
HKomatsu et al.ll2009l ). 
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Figure 7. Top three eigenvectors (EVs) chosen from the nor- 
malized covariance matrix in the Mr < —21.5 galaxy sample. The 
sign of the EV is arbitrary. The first EV (solid black) shows equal 
weights for all bins. The second (dashed red) and third (dotted 
blue) EV display the configuration difference between perpendicu- 
lar and co-linear triangles as well as the scale variation as the scale 
of the third side increases. 
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Figure 8. Like Figure[7]but for the -21.5 < Mr < -20.5 galaxy 
sample. The top three eigenvectors (EVs) chosen from the nor- 
malized covariance matrix. The sign of the EV is arbitrary. The 
first EV (solid black) shows equal weights for all bins. The second 
(dashed red) and third (dotted blue) EV display the configuration 
difference between perpendicular and co-linear triangles as well as 
the scale variation as the scale of the third side increases. 



we are most confident with using the local bias model). 
Given our analysis of the relative bias, we use the larger 
B values where we constrain C = 0. We present these 
estimates in Table HI 

5. EIGENVECTORS OF THE 3PCF COVARIANCE 
MATRIX 

A point that is often overlooked is that the covariance 
matrix itself is a measurement of clustering rather than 
simply a means of quantifying uncertainty. It exhibits 
increased sensitivity to higher order terms (for a concise 
review see Szapudi 2009) with the covariance of the 2PCF 
being leading order sensitive up to fourth order, and the 
3PCF up to sixth order. 

We investigate the structure of the normalized covari- 
ance matrices by examining the eigenvectors (EVs), or 
principal components, obtained by a singular value de- 
composition. The EVs are contained in the U and V 
matrices from ([TSl) and (fT6|). The first EV is associated 
with the largest singular value (SV) , and accounts for the 
largest variance in the normalized covariance matrix (i.e. 
most of the observed structure); the second EV is the 
next largest SV and so on. If the covariance matrix re- 
solves predominately "true" signal, the first EVs should 
characterize this structure whereas the lower ranked EVs 
encapsulate noise. While the amplitude of the EVs are 
not significant without the corresponding SV, they do 
represent the variation between bins in an orthogonal 
basis where the full covariance is a simple linear combi- 
nation. 

We show the top three EVs for the BRIGHT and 
LSTAR galaxy samples in Figures [7] and [51 respectively. 
In all cases there appear to be consistent features in 
the eigenmodes. The first EV represents weighting all 
bins equally. Typically, the second EV highlights the 



difference between "perpendicular" and "co-linear" con- 
figurations. Finally, the third EV tracks a roughly 
monotonic change from small to large 0, possibly ac- 
counting for the scale difference of the continually in- 
creasing third side of the triangle. We point out that 
this structure evident in observational galaxy samples 
agrees well with theoretical predictions from simulations 
in Gaztafiaga & Scoccimarro (2005). Remember, these 
EVs are obtained by deconstructing just the normalized 
covariance matrix. 

We do not always see a clear separation between the 
second and third EVs. As the full structure is a linear 
combination of all modes, the configuration dependence 
and scale variation effects could be combined. The SVs of 
the two effects are essentially equivalent for our measure- 
ments, making their numerical distinction in the SVD 
somewhat arbitrary. This is not a concern, as it appears 
that the linear combination of these two EVs is consis- 
tent with our interpretation (even if they are mixed). 
Less significant EVs show less coherent structure, con- 
sistent with noisy modes in the covariance (as we would 
expect). 

By examining the EVs of the covariance matrices, we 
note structure consistent with measurements of the re- 
duced 3PCF (see Figures [3] and HJ. Observing this struc- 
ture provides supporting evidence that we have signal 
dominated estimates of the covariance matrix. This jus- 
tifies our approach of using a combination of the most 
significant eigenmodes in a quantitative comparison to 
galaxy-mass bias models, as we did in ^ 

6. QUALITY OF ERROR ESTIMATION 

We rely heavily on the structure of the error covariance 
matrix for constraints on galaxy-mass bias. We noticed 
the observed structure in the covariance is qualitatively 



Redshift r^.r^ = 9:18 Mpc 



Projected r^r^ = 9:18 Mpc 
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Figure 9. We compare the Icr absolute (diagonal) errors of the reduced 3PCF obtained by using different methods of estimation: 
independent mock catalogs or jackknife resampling as denoted. These measurements correspond to the BRIGHT galaxy sample. 



similar to clustering measurements (in fj5|) , but it remains 
unclear if this structure is affected by our jackknife re- 
sampling estimation of the covariance. Higher orders add 
complexity and increased sensitivity to systematics, even 
with a "ratio statistic" such as the reduced 3PCF where 
the error sensitivity is canceled to first order. We must 
investigate the error resolution of jackknife resampling on 
the 3PCF, as tests on the 2PC F in angular correlation 
functions (iScranton et al.l 2002 ) or redshift space SDSS 



(|Zehavi et al.ll2002L I2005r 



2010 ) can not be assumed to 



be sufficient. 

An alternative method to estimate errors uses a se- 
ries of independent realizations of artificial galaxies, ide- 
ally created to match observational limitations such as 
the volume and geometry of the SDSS galaxy samples. 
We created 49 independent galaxy mock catalogs based 
on independent A'^-body simulations that have appro- 
priate resolution to match the BRIGHT SDSS sample 
{Mr < —21.5). We use these independent mocks to 
estimate errors and compare with those obtained from 
jackknife re-sampling of the data. This exercise should 
provide an idea how effective jackknife re-sampling is for 
resolving the errors on the 3PCF. The BRIGHT galaxy 
sample has the lowest number density with the least 
number of galaxies over the largest volume. To help pro- 
tect against undersampled measurements due to low bin 
counts (see discussion in IMcBride etlll[2010h . we re- 
strict the comparison to the configuration dependence of 
the larger triangles (ri = 9 h^^Mpc sides). 

We estimate the covariance matrix for the BRIGHT 
sample using different numbers of jackknife samples on 
the data, specifically using 15, 30, 49 and 105 jackknife 
regions. Again, these jackknife sample are created from 
the observational data and not from the mock galaxy 
catalogs, where each jackknife region is selected to main- 
tain equal unmaske d area (same method as detailed in 
IMcBride et alllMnh . Since we measure 15 bins for Q{0), 
we require at least 15 jackknife samples to prevent a 
singular covariance matrix. We use twice this number 



(30) and then use the same number as the number of 
mocks (49). The final value corresponds to the number 
of unique elements in the symmetric covariance matrix: 
15(15 — l)/2= 105. We caution that as we increase the 
number of jackknife samples, we decrease the respective 
volume of each jackknife region which might subtly bias 
jackknife estimates (e.g. underestimate cosmic variance). 

First, we investigate the magnitude of the absolute (di- 
agonal) errors of the reduced 3PCF. Since we use nor- 
malized covariance matrices, differences in the absolute 
errors might not be noticeable in the covariance struc- 
ture. The 1(7 absolute errors are shown in Figure |9l We 
see little difference between any of the methods, and the 
uncertainty typically ranges between 0.1 and 0.15. 

For each of our methods, we estimate the normalized 
covariance in both redshift and projected space, as de- 
picted in Figure [TOl and include the distribution of resid- 
uals. We note that jackknife re-sampling methods ap- 
pear to underestimate the correlation in all cases, but 
the general structure looks comparable. More samples 
generally produce a smoother, and more correlated, co- 
variance matrix. However, not even the 105 jackknife 
sample estimate reproduces the correlation in 49 mocks. 
We consider the distribution of residuals an important 
metric in evaluating the reliability of the resulting co- 
variance, which we include in Figure flOl Ideally, the co- 
variance matrix accounts for all "connections" between 
bins only if the residuals are reasonably Gaussian. We 
notice a skew in several of the jackknife re-samplings, 
with a tail ext e nding to lower values. As discussed in 
IMcBride etall (|2010t i. this is a consequence of cosmic 
variance within jackknife samples. A few rare structures 
affect the 3PCF; when they are excluded by an a jack- 
knife region the Q{0) of the entire sample drops. The 
mock estimate shows a slight skew in the positive direc- 
tion from the same effect. In mocks, when a rare struc- 
ture exists in the probed volume then the 3PCF rises 
producing a rare high measurement. 

The eigenmode analysis we utilize relies on signal being 
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Figure 10. We present the normalized covariance matrices and residuals of the error estimation for large triangles (9 — 27 /i~^Mpc). The 
left and right columns pertain to redshift and projected space respectively. We estimate errors using 15, 30, 49, and 105 jackknife regions 
and compare with results from 49 mocks from independent A^-body simulations. The solid, dashed and dotted contours in the normalized 
covariance correspond to values of 0.70, 0.85, and 0.99, respectively. 
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Figure 11. Top three eigenvectors (EVs) chosen from the normal- 
ized covariance matrix for the different error estimates for Qz{0) 
and Qproji^)- The sign of the EV is arbitrary. The first EV (left 
panels) shows equal weights between all bins. The second (middle) 
and third (right) EV display the configuration difference between 
perpendicular and co-linear triangles, as well as the scale variation 
as the scale of the third side increases. 
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Figure 12. The singular values (SV), or eigenvalues, obtained 
from the singular value decomposition (SVD) of the normalized 
covariance matrix for each of our error estimates. Larger values of 
the SV correspond to more statistically significant eigenmodes in 
the structure of the covariance. 



the dominant contribution to the structure of the covari- 
ance matrix (as opposed to noise). Noise is commonly 
expected to be an independent or diagonal contribution. 
Similar to iJSl we examine the eigenvectors (EVs) of the 
covariance matrix to provide insight into the structure. 
By using the singular value decomposition (SVD), the 
eigenvectors are ordered by largest to least amount of 
variance explained in the covariance matrix. 

The first three EVs are shown in Figure [11] for both 
redshift and projected space. Similar structure appears 





Figure 13. The signal-to-noise ratio for each eigenmode, ordered 
in terms of importance. The total signal-to-noise of a measurement 
is calculated by adding each individual eigenmode in quadrature. 
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Figure 14. The cumulative signal-to-noise ratio for each eigen- 
mode, ordered in terms of importance. The cumulative total is 
calculated by summing in quadrature the more significant modes. 



in each of them, which we interpret as follows. The 
first EV represents the general measurement, with all 
eigenmodes equally weighted. The second EV shows 
the difference between "collapsed" and "perpendicular" 
configurations. Finally, the third EV represents a scale 
dependence as the third side of the triplet ranges be- 
tween 9 /i-^Mpc at 6' to 27 /i"^Mpc at 6* ~ tt. 
In some of the estimates, the shapes of the second and 
third EVs appear either combined or transposed. Since 
the full measurement is a linear combination of all EVs, 
this lack of separation makes sense. In these cases, 
the statistical significance the two EVs remain similar. 
T his interpretation of the structure follows the analysis 
bv lGaztanaga fc Sco ccimarrol (|2005[ ) for V-body simula- 
tions. The less significant eigenvectors (which we do not 
show) appear random, with the lowest being contribu- 
tions from noise or numerical instabilities. We identify 
the significance of the eigenmodes by inspecting the sin- 
gular values (SVs) shown in Figure [121 The SV can be 
understood as an "importance weighting" of each eigen- 
mode, and the figure shows a rapid decline of significance 
for each eigenmode. The first three eigenvectors cumu- 
latively account for over 99.9% of the variance in the 
normalized covariance matrix. 

The signal-to-noise ratio (S/N) of each eigenmode is 
shown in Figure [T51 as calculated by (HZ]). The mocks 
in both redshift and projected space depict a slow de- 
cline in S/N over the first few eigenmodes, supportive or 
our interpretation of relative significance. This trend is 
not as clear in the jaekknife estimates for redshift space, 
although it appears consistent in projected space. We 
see the first half of the modes appear resolved, with 
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Figure 15. We compute the compatibility of the subspace contained in a series of eigenvectors sucii tiiat the y-axis can be interpreted 
as the fractional "match" between the spaces the eigenvectors probe. A value of 1.0 means no mismatch in the space they probe, and 0.0 
means no overlap (i.e. orthogonal). The left panel uses the eigenvectors of the redshift space covariance matrix determined by the 49 mocks 
as the reference value. The right plot does the same as the left, but uses measurements in projected space. The comparison is cumulative 
(eigenmode 3 means the sum of the first 3 modes). 
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Figure 16. We use the subspace comparison of eigenvectors to 
estimate the difference in space probed between similar numbers of 
eigenmodcs in redshift and projected covariance matrices for each 
of the error estimates. 



well behaved S/N. For the least significant eigenmodes, 
the noisier error estimates using fewest jackknife sam- 
ples show unrealistically high S/N ratios (especially in 
the case of 15 jackknife regions). The total S/N would 
increase dramatically and artificially if we included these 
noise dominated modes. In these cases, using the full 
covariance (i.e. including all modes) would be a mis- 
take. To make the point clearer, we examine the cumu- 



lative S/N ratio in Figure [14] where we identify rapid 
upturns in the total S/N as an artificial consequence of 
noise. Several curves in Figure [M] do not appear prob- 
lematic with this test, and show steady behavior across 
all modes. The amplitude of the S/N ratio between 49 
mocks and the 105 jackknife samples show consistency, 
but the S/N ratio does not appear to be a monotonic 
change with the number of jackknife samples which sug- 
gests a complex relationship between the best S/N and 
an optimal number of jackknife regions. 

We can compare the subspace that a set of eigenvec- 
tors probe between two err or estimat e s. Th e formalism 
is the same as discussed in lYip et al.l (|2004|, see section 
4), which results in a fractional "compatibility" between 
a collection of eigenvectors. Intuitively, this is the matrix 
equivalent of the vector dot product, where two orthog- 
onal unit vectors would have a vector subspace of (no 
compatibility) and two identical unit vectors would result 
in 1. We use the covariance of the 49 mocks as "truth", 
and test the fractional compatibility of the jackknife esti- 
mates for covariance in Qz{0) and Qprojifl) shown in Fig- 
ure ll5l When all the eigenmodes are considered, the sub- 
space becomes the full space and the comparison yields 
unity by construction. We notice the projected measure- 
ments never appear more discrepant than 75%. After 
the first few eigenmodes, redshift space shows a similar 
agreement. With the exception of the 15 jackknife sam- 
ple estimate, the 3 eigenmode mark appears 90% com- 
patible or better in all cases. This quantifies our argu- 
ment of the top three EVs in Figure [TTl where the second 
and third eigenvectors appear different (predominantly 
in redshift space), but their linear combination remains 
consistent with each other. Remember, this comparison 
only considers the compatibility of the direction of each 
eigenvector, and not their relative strengths (i.e. SVs). 
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We evaluate the subspace compatibility on the normal- 
ized covariance matrix between redshift and projected 
space estimates. For each method, we show the fractional 
comparison in Figure 1161 The mocks estimates show the 
most compatibility across all eigenmodes, where showing 
agreement at ^ 85% or better. With all estimates, we 
find that the combination of the first three eigenmodes 
remains a compatible subspace, above 90%, if we again 
exempt the 15-jackknife sample estimate (which shows 
less than 70% compatibility). 

We caution that the resolution of errors and the choice 
of binning scheme relate in a non- trivial manner, which i s 
discussed in additional detail bv iMcBride et al.l ()2010D . 
We chose "large" bins (fiducial scheme with / = 0.25) 
to ensure a smooth, signal dominant structure in the co- 
variance matrix. Overall, this error comparison supports 
our claim that accurate results can still be obtained even 
with less-than-optimal error estimation such as jackknife 
re-sampling. 

7. DISCUSSION 

We utilize the the configuration dependence of the 
3PCF in redshift and projected space to constrain 
galaxy-mass bias parameters in the local bias model. 
We find that galaxies are biased tracers of mass, with 
brighter galaxies corresponding to increased bias. These 
results are consistent with deta iled analysis of SDSS 
galaxies from the 2PCF (jZehavi ct al. 2005, 2010) which 
quantifies how bias increases clustering for brighter 
galaxy samples. Our results indicate that a linear bias 
model yields reasonable approximations to t he observa- 
tions, in agreement with iHikage et al.l ()2005[ ). However, 
a non-linear bias model produces slightly better agree- 
ment, and yields lower reduced chi-square values {xt i^i 
Tables [5] and [3]) . We notice a strong correlation between 
linear and quadratic bias, as expected from inspection of 
pT|) . and consistent with measurements of SPSS galax- 
ies using the bispectrum (|Nishimichi et al.ll2007D . We 
find that our redshift space measurements predict signifi- 
cantly negative quadratic bias with a linear bias near one. 
This effect was se en in a similar an alysis conducted on 
2dFGRS galaxies (jGaztafiaga et al.l[2005, '). Interestingly, 
we find projected measurements suggest a larger linear 
bias with near zero quadratic bias for the same samples, 
suggesting a possible systematic effect from redshift dis- 
tortions in this simple bias model. 

We examined the relative bias in §4.31 We find sup- 
porting evidence that the brighter galaxy sample is 
a more biased realization using both the 2PCF and 
3PCF, consistent with other analyses of SDSS data 
(jZehavi et al.l 120021 120051 12010D . Relative bias provides 
a consistency check on the "absolute" galaxy-mass bias 
parameters we constrain, suggesting a combination of lin- 
ear and quadratic bias terms are consistent with obser- 
vations. However, the relative bias of the 2PCF suggests 
that our two parameter bias model fits underpredict the 
value of linear bias necessary to explain the observations. 
Again, we see a hint that constraints from projected mea- 
surements appear to be less affected - although we cau- 
tion that this trend has weak statistical significance given 
the larger uncertainties in projected space. 

We obtain reasonable projections for erg by using our 
linear bias values from fits on Q{9) in conjunction with 
the 2PCF. We estimate the values of as to be between 



0.83 and 1.13 based on the BRIGHT {Mr < -21.5) and 
LSTAR (-21.5 < Mr < -20.5) galaxy samples. The 
values we obtain are contingent on a specific model of 
mass clustering, where we have chosen to use A'^-body 
simulations (specifically the Hubble Volume ACDM re- 
sults), and redshift distortions (which we include through 
velocity information to distort particle positions in the 
HV simulation). For comparison, constraints of erg from 
a joint analysis of the cosmic microwave background 
(CMB), supernova data (SN) and baryon acoustic os- 
ciUations (BAO) find as = 0.82 (Komatsu et al. 2009). 
Our lower values are in good agreement with these con- 
straints. Our high end values appear too large, but 
our results are in reasonable agreement with an anal- 
ysis of a related statistic, the monopole moment of the 
3PCF, where they find best fit as values between 0.9 and 
1.07 (see Table 3 in 'Pan & Szapudi"2005) using 2dFGRS 
galaxies (CoUess et a l. 2001). Although the value of erg 
we obtain is comparable with results from 2dFGRS, the 
specific bias values will not be, as the 2dF targets a difi'er- 
ent galaxy selection than our SDSS samples. If we under- 
estimate the value of linear bias, effectively B here, (|23|) 
shows that the implied value of erg will be overestimated. 
This might explain the larger values of our estimates in 
comparison to WMAP analyses. Our projections for erg 
use clustering measurements between 9 — 27 h^^Mpc 
and exploit only the configuration dependence of Q{9). 
This is a much smaller slice of data that is significantly 
different than either the monopole measurement (which 
utilizes a larger range of scales without configuration de- 
pendence) or WMAP results (that combines a immense 
amount of data from both CMB and LSS analyses). We 
do not intend this analysis to complete with these con- 
straints, but rather to help illuminate the role of galaxy- 
mass bias in future constraints of cg using the 3PCF. 

Understanding the properties of measurement errors 
and the impact of empirical methods of estimating the 
covariance is a critical component necessary for quan- 
titative constraints. Recent results have done compar- 
isons on lowe r orde r statistics, such as the work by 
iNorberg et al.l (|2009f ). We compared several properties 
of 3PCF covariance matrices estimated from jackknife 
re-sampling to those constructed from many realizations 
of independent galaxy mock catalogs. While we noted 
some concerning discrepancies, we found these typically 
affected only the least significant eigenmodes. We found 
many similarities between the covariance estimates, in- 
cluding physical descriptions for the first three eigen- 
modes which account for an overwhelming majority of 
the variance. We established the need to trim noisy, un- 
resolved modes from the covariance. When trimmed, and 
the eigenmode analysis is properly utilized, we noted only 
a few significant differences, mostly in the case of 15 jack- 
knife samples. We conclude that our use of 30 jackknife 
samples does not significantly affect our analysis. 

8. SUMMARY 

We analyze measurements of the configuration depen- 
dence of reduced 3PCF for two SDSS galaxy sam ples 
that were first presented in IMcBride et al.l (|2010l ). In 
both redshift and projected space, we characterize the 
galaxy clustering differences with those predicted by the 
non-linear mass evolution in the ACDM Hubble Volume 
simulation. Here, we summarize our main results: 



We demonstrate that brighter galaxies remain a 
more biased tracer of the mass field by constraining 
the linear and quadratic galaxy-mass bias parame- 
ters using a maximum likelihood analysis on scales 
between 6 and 27 h~^Mpc. Conservatively using 
scales above 9 h~^Mpc, the BRIGHT sample is 
biased at greater than 2a and the fainter LSTAR 
shows no significant bias, in generally agreement 
with exp ectations frorn prev i ous a nalyses of SDSS 
galaxies (jZehavi et all I2005L l2010f ). The bias pa- 
rameters and their significance are summarized in 
Table H 

We resolve the degeneracy between the linear and 
quadratic bias terms, which helps to explain the 
weak luminosity dependence observed in the re- 
duced 3PCF. 

We find a linear bias model appears sufficient to ex- 
plain the measurements of the 3PCF by re-fitting 
the linear bias while constraining the quadratic bias 
at zero (results reported in Table [3]). However, we 
find the two parameter fit is preferred in our likeli- 
hood analysis, as it yields a lower chi-square in the 
best fit value. 

The relative bias between samples of different lumi- 
nosities (which is independent of the mass predic- 
tions), as well as the cosmological implications for 
values of erg , show general consistency with previ- 
ous analyses. Inspection of our results suggest that 
the linear bias values obtained without a quadratic 
bias term are preferred. This suggests that two- 
parameter bias constraints might underpredict the 
linear bias. 

We decompose the structure of the normalized co- 
variance matrix as an alternative view into clus- 
tering properties of our samples. The eigenvec- 
tors of the first three dominant modes show co- 
herent structure consistent with variations seen in 
the Q{0) measurements, supporting our claim that 
the covariance is signal dominated and sufficiently 
resolved. 

We find that jackknife re-sampling methods cannot 
reproduce the correlation seen in the a 3PCF co- 
variance matrix estimated from many realizations 
of mock galaxy catalogs. By performing a detailed 
comparison of the properties and structure of the 
errors, we identify that noisy, unresolved modes in- 
troduce significant discrepancies. We find that us- 
ing an eigenmode analysis can mitigate the differ- 
ences and conclude that our analysis should not be 
significantly affected by less-than-ideal methods of 
error estimation. 

Comparing results between redshift space and pro- 
jected measurements implies a potential systematic 
bias on values from the redshift space analysis when 
scales below 9 h^-^Mpc are included, which have 
been utilized in other comparable analyses. Since 
the small scale measurements contain more con- 
straining power than larger scales, they drive the 
likelihood analysis even when larger scales are con- 
sidered. 



• On scales above 9 /i^^Mpc, the statistical signif- 
icance of constraints from redshift space analy- 
ses appear stronger than those found in analyses 
of projected measurements. We attribute this re- 
sult to the increased uncertainties of the projected 
3PCF, which mixes in larger scales (with larger er- 
rors) due to the line-of-sight projection. When con - 
sidered with the results of iMcBride et al.l (|2010D . 
which finds the projected 3PCF recovers configu- 
ration dependence at small scales lost in redshift 
space, a combination of redshift space analysis at 
large scales and projected measurements at small 
scales would form a nice complement in future anal- 
yses. 
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Figure Al. Analogous to the Galaxy-Mass bias constraints in ijjl we show the constraints on B and C using the same data for our 
fiducial binning scheme with fractional bin-size of / = 0.1 (solid red contours) and larger / = 0.25 (block dotted). On the left, we neglect 
any overlap as well as the the covariance and assume independent diagonal errors while and use the full 15 bins. On the right, wc utilize 
the full covariance and only fit the dominant modes in an eigenmode analysis. The contours correspond to the Itr and 2a confidence levels 
from the Ax^ surface. We use the LSTAR galaxy sample. 
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APPENDIX 
A. EFFECTS OF BINNING 

As an example of the effects of bin-size on galaxy-mass bias constraints, we re-analyze the LSTAR galaxy sample 
(—21.5 < Mr < —20.5) using the two fractional bin-widths: / — 0.1 and / = 0.25. First, we ignore the structure 
of the covariance matrix and show constraints using all the bins assuming perfect independence shown in left panel 
of Figure lAll While unphysical, this illustration allows one to probe the effect of shape differences in the 3PCF 
measurements without considering the resolution of the covariance matrix. Since larger bins smooth the configuration 
dependence, we expect a larger degeneracy between B and C, which is apparent. We see the best fit values (symbols) 
stay within the respective la contours, but just barely. Remember, this approach uses all modes and the exact same 
input data, suggesting that binning can result in a la systematic bias. 

For the right panel of Figure IA1[ we consider the full covariance as well as improvements obtained by using the 
eigenmode analysis in the galaxy-mass bias constraints. First, we notice the error contours appear less stretched, 
in accord with our expectations of using a non-diagonal covariance matrix. In most cases (excepting for ri = 
6 h~^Mpc), the area of the contours appear of equal size or even decreased for the larger / = 0.25 measurements in 
contrast to the diagonal case. This makes sense, as the lower variance measurements of / = 0.25 appear better resolved 
as long as there are enough remaining modes to constrain two parameters. The best fit values appear discrepant, 
especially at the lower scales (6 — 18 h~^Mpc) where they disagree at more than a la significance. While this causes 
some concern, it is not as drastic as the diagonal case. As the eigenmode analysis trims modes, it excludes information 
and the same input data produces a different statistical representation. In light of this effect, a la difference becomes 
a statistical difference of analysis rather than a significant systematic effect. 

In summary, we find lower galaxy-mass bias parameters with larger bin-widths, a potential artificial bias on the 
galaxy-mass parameters due to over-smoothing. Since we gain very little additional constraining power with the 
/ = 0.25 bin-width, we argue the / = 0.10 bin-width represents the more conservative choice. Although the / — 0.10 
scheme represents smaller bins, they are still quite large and adequately resolve structure in the covariance. 
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